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Comparison between a quantum kinetic theory of spin transfer dynamics in Mn doped 
bulk semiconductors and its Markov limit for non-zero Mn magnetization 
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We investigate the transfer between carrier and Mn spins due to the s-d-exchange interaction in 
a Mn doped bulk semiconductor within a microscopic quantum kinetic theory. We demonstrate 
that the spin transfer dynamics is qualitatively different for components of the carrier spin parallel 
and perpendicular to the Mn magnetization. From our quantum kinetic equations we have worked 
out the corresponding Markov limit which is equivalent to rate equations based on Fermi’s golden 
rule. The resulting equations resemble the widely used Landau-Lifshitz-Gilbert-equations, but also 
describe genuine spin transfer due to quantum corrections. Although it is known that the Markovian 
rate description works well for bulk systems when the initial Mn magnetization is zero, we find large 
qualitative deviations from the full quantum kinetic theory for finite initial Mn magnetizations. 

These deviations mainly reflect corrections of higher than leading order in the interaction which are 
not accounted for in golden rule-type rates. 
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Diluted magnetic semiconductors (DMS) have been 
studied intensively in the past decades, since they com¬ 
bine the versatility of semiconductors with the spin de¬ 
gree of freedom, which promises future applications in 
spintronics^^®. The magnetic properties of DMS arise 
from the s/p-d exchange interaction^’®’^ between carriers 
and magnetic impurities, which typically consist of Mn 
ions acting as localized spin | systems. Especially for 
short timescales and high Mn doping concentrations the 
exchange interaction can dominate the spin dynamics®’®. 
The description of the resulting spin transfer dynamics 
in DMS is usually based on rate equations, where the 
rates are computed using Fermi’s golden rule®’^®. The 
standard derivation of the golden rule involves a Markov 
approximation®’^^ and is perturbative with respect to the 
exchange coupling constant. In Ref. 12 a projection oper¬ 
ator method was applied to derive spin relaxation rates 
for DMS quantum wells. There, also a Markovian as¬ 
sumption as well as a perturbative argument were used. 
Another approach to the description of the macroscopic 
magnetization dynamics is the use of the phenomenolog¬ 
ical Landau-Lifshitz-Gilbert equations^®’^^. 

Recently, starting from a Kondo-like interaction 
Hamiltonian a density matrix approach based on correla¬ 
tion expansion was developed®® in order to describe the 
spin dynamics in the ultrafast regime. Until now, this 
quantum kinetic theory (QKT) has only been applied 
to the case of an initially zero Mn spin. There, it has 
been found that in three dimensional systems, the time 
evolution of the carrier spin is exponentially decreasing, 
where the decay rate coincides with its value according to 
Fermi’s golden rule®®. The latter was shown by perform¬ 
ing the Markov limit (ML) of the QKT using only terms 
in second order of Jsd- In lower dimensional systems, 
excitation conditions can be found, where significant dif¬ 
ferences between the ML and the QKT become visible 
although the memory induced by the exchange interac¬ 
tion is orders of magnitude shorter than the timescale 
for the evolution of the carrier and Mn dynamics®®. In 


particular, quantum kinetic effects are most pronounced 
when suitably tuned oppositely circular polarized two- 
color laser pulses are used for the excitation®®’. 

In this article, we study the spin dynamics of conduc¬ 
tion band electrons in a bulk ZnMnSe semiconductor for 
the case of a non-zero initial Mn spin where electron spins 
can process around the Mn magnetization. It turns out 
that the spin transfer dynamics that is superimposed to 
the precession is qualitatively different for electron spins 
aligned parallel or perpendicular to the Mn magnetiza¬ 
tion. Starting from our quantum kinetic equations we 
derive the corresponding Markov limit for finite Mn mag¬ 
netization. The resulting equations can be interpreted 
as modified Landau-Lifshitz-Gilbert equations. Assum¬ 
ing Mn concentrations much larger than the itinerant 
electron density analytical solutions of these Markovian 
equations are presented. The resulting analytical expres¬ 
sions also exhibit a different dynamics for perpendicular 
and parallel spin transfer, which, however, quantitatively 
and qualitatively disagrees with the prediction of the full 
QKT. Here, the failure of the Markovian approach can 
be traced back to contributions of higher than leading 
order in the exchange coupling constant. 

Outline of the paper; In a first step, we briefly summa¬ 
rize the QKT®® that was used as a basis for our numerical 
calculation and introduce the model used in this paper. 
Then, we derive the Markov limit of the QKT along the 
lines described in Ref. 16 for an initially zero Mn magne¬ 
tization (S), but allow for a finite value of (S) and an ar¬ 
bitrary angle between the conduction band electron spin 
and the Mn spin. In a subsequent section we present nu¬ 
merical results of our QKT for the spin transfer dynamics 
of the parallel and perpendicular components and com¬ 
pare them with the ML. The analytical solution of the 
ML equations in combination with a rearrangement of 
the contributions to our QKT allows for a clear physical 
interpretation of the pertinent source terms. By selec¬ 
tively studying the impact of different source terms we 
are able to demonstrate the importance of contributions 
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of higher than leading order in the coupling constant. 


I. QUANTUM KINETIC EQUATIONS 

In Ref. 15, a quantum kinetic density matrix approach 
for the spin dynamics in Mn doped semiconductors was 
developed starting from the Hamiltonian: 


H — Hq + Hsd + Hpd + Hern, ( 1 ) 

where Hq describes the single particle band energies, Hgd 
accounts for the exchange interaction between the s-type 
conduction band electrons and the spins of the d-type 
electrons of the Mn dopands while Hpd stands for the 
interaction of the latter with p-type holes. Finally, idem 
comprises the dipole coupling to an external laser field. 
The exchange interactions Hgd + Hpd as well as the ran¬ 
dom spatial distribution of Mn atoms give rise to a hi¬ 
erarchy of higher order correlation functions. In order 
to obtain a finite set of dynamical variables a specially 
adapted correlation expansion has been worked out in 
Ref. 15. 

Since the aim of the present paper is to investigate the 
spin transfer between conduction band electrons and Mn 
dopands, the model can be reduced to: 

H = Ho + H,d- (2) 

Hq now accounts only for electrons in a single spin de¬ 
generate conduction band: 

d^o = (3) 

Ik 

where cjj^ (cik) are the creation (annihilation) operators 
of conduction band electrons with k-vektor k and spin 
index I = 1,2. For simplicity we shall assume parabolic 

jr 2 

bands E\^ = with an effective mass m*. The ex- 

I 


change interaction is given by^®’^®: 

Hsd = JsdY.^I■mr^-'RI), (4) 

li 


where Jsd is the exchange constant and S/ (s®) are op¬ 
erators for the spin of the Mn atom (conduction band 
electron) in units of h at the position R/ (r^). As in 
Ref. 15 we assume an on average spatially homogeneous 
distribution of Mn positions R/. 

According to the analysis in Ref. 15 the relevant dy¬ 
namical variables for this reduced model are: 

= (4iki'^*2ki), (5a) 

MT, = {PLn,), (5b) 

(5c) 

(5d) 

where := |d,ni)(d, 712 ! describes the spin state of 

the Tth Mn ion (n = —1,...,|). The expectation value 
represented by the brackets involves a quantum mechan¬ 
ical average as well as the disorder average over the ran¬ 
domly distributed Mn positions. and are the 

electron and Mn density matrices. and are 

the correlated parts of the corresponding density matri¬ 
ces, i.e., in these quantities all parts that can be factor¬ 
ized into products of lower order correlations functions 
are subtracted from the expectation values. The explicit 
but lengthy definitions of be found 

in Ref. 15. 

It turns out that the resulting equations of motion can 
be simplified by introducing the following new correlation 
functions: 


^I2n2k2 

^Ziniki 






( 6 ) 


Rewriting the equations of motion from Ref. 15 in terms 
of these functions we obtain: 


^ ^ + i ^ (Snn.QLk*^' - ^n.nQtX) 


V 


k nil' 


V ■ 


= Jmiu ^ MX{su,Cjl^ - + -^ E - ^hiQ\:X) 

nn'l ^ k 

+ Ak, - + b\l2kf + b\TSf'^ 


with source terms 


(7a) 

(7b) 

(7c) 


Z2n2k2 ^ 
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^Zk2" 


(S„„,M"= - Sn,nM:^ 


), (7d) 
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nil' k 


nn'l 
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L^2^2k2 


(7e) 


f,hn2k.2^^^ —J / _L [o O,, /Q*2"2k2 _ C „ ,/Q*"k ] 


, y 

nl ^ k 


,12^2^2 


- J.d E E [^nnMT'^^ - ^n2nQtt‘] + ^ ^ ^'z! [Snn. " S^.nQf.tkJ |, 

nil' ^ k k 


J-2^2^2 ^^^-2 


(7f) 


where S„j ^„2 are the Mn and electron spin matri¬ 

ces, V is the Volume of the DMS and riMn = is the 
density of the Mn ions. We have subdivided the sources 
on the r.h.s. of Eq. 7c for later reference. The physical 
meaning of these terms and their respective importance 
will be discussed later. 

In order to study the dynamics of the spin transfer 
we consider initial conditions where the electrons are ini¬ 
tially spin polarized and the Mn magnetization corre¬ 
sponds to a thermal distribution while the correlations 
Qhnlkl assumed to be zero. This is a situation typ¬ 
ical for a system, immediately after an ultrafast optical 
excitation has induced a finite electron spin polarization. 


II. MARKOV LIMIT 


It turns out to be instructive to derive the Markov limit 
of our QKT, first of all, because this greatly simplifies the 
theory as the higher order correlation functions are for¬ 
mally eliminated in favor of the variables of most inter¬ 
est, i.e., the electronic densities and spins. Furthermore, 
the Markov limit provides a relevant reference for our 
QKT. In particular for bulk systems it has been found 
previously^® that the memory of the exchange interac¬ 
tion is short and therefore it is tempting to think that 
the Markovian equations should yield valid results in our 
case. 

In order to be able to work out the Markov limit start¬ 
ing from Eqs. (7), we follow the procedure that in Ref. 16 
led to rates in accordance with Fermi’s golden rule and 
neglect in a first step the source terms of higher than 
leading order in the exchange coupling Jgd- Due to the 
initial condition = 0 the correlations are 

of first order in Jgd and thus we see from Eqs. (7) that 


and of second order in J^d and 

yield third order contributions to the electron spin dy¬ 
namics. Thus, we keep in Eq. (7c) only the first order 

term . This allows us to formally integrate the 

correlations: 


t 

Q\TnX\it) = ^ /(8) 

0 

with frequency Wk = -^ = 5 ^. Substituting Eq. ( 8 ) 
back into the equations for and we have to 

perform a fc-summation which due to interference result¬ 
ing from the fe-dependent phases -t) 

a finite memory. The Markov limit is established by as¬ 
suming that the sources change on a much slower 

timescale than the memory and can therefore be drawn 
out of the integral. The memory has been found to de¬ 
cay on a fs-timescale while the spin dynamics evolves on 
a ps-timescale^®. Therefore, the lower limit of the inte¬ 
gral can be extended to —00 resulting in the following 
approximation for the correlations: 


0 

—00 


_ lkZ2n2k2^C^'| 
- fc^Zmiki VI 



Wki) -75 - 

Wk2 


i 

-a;ki 


(9) 


where V denotes the Cauchy principal value. 

Starting from Eq. (7b) for the electron density C'jjki 
can set up an equation of motion for the average electron 
spin (ski) = ®q/ 2 ^Ziki state with k-vector 

ki. Feeding back the correlations from Eq. (9) 

into these equations we finally obtain: 
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d , , J sd'^Mn / jr^\ / , J sd^Mt] 




+ 7r(5(wki - “^k) 


(S) 


4{ski)^ - + 2nki 


+ ((sk) X ((ski) X (S))) + 


(S X (S X (sk))) + (((sk) X S) X S) 


( 10 ) 


Applying the same procedure to the electron occupations 
riki = X)z ^zki ^ given k-vector ki we find that on 
this level of theory Uki is time independent. It should be 
noted, that in the full QKT this is not the case. Instead it 
was shown in Refs. 15-17 that redistributions in fc-space 
take place which are responsible for a number of features 
of the magnetization dynamics that are not expected in 
the Markovian theory. 

The different terms in equation (10) can easily be in¬ 
terpreted. The first term describes the precession of the 
electron spin in an effective magnetic field due to the Mn 
magnetization (S), which is also the result of a mean- 
held calculation^®. The second term represents a renor¬ 
malization of the precession frequency that depends on 
the density of states and therefore on the dimensionality 
of the system as well as the k-vektor, which can possibly 
lead to dephasing of the electron spin. 

The magnitude of the renormalization for a bulk semi¬ 
conductor can be estimated in the continuum limit by 
approximating the Brillouin zone (BZ) as a sphere with 
radius ksz and assuming a parabolic band structure as 
follows: 

* 2 

0 

^ V 

where |(s^| jg the mean-held precession fre¬ 

quency. The order of magnitude of the integral on the 
r.h.s. of Eq. (11) can be determined by noting that the 
optically excited carriers occupy only very few states near 
the center of the BZ and therefore for the most part of 
the BZ rik ~ 0 holds which also implies ^ ^ 1 for the 
occupied states. Approximating rik ~ 0 and « 0 the 
integral yields the value kBz- For the parameters used 
in our study (see below) the renormalization is estimated 
in this way to be of the order of ~ 1% of the mean-held 
precession frequency^'^. 

The third term in Eq. (10), which is proportional to the 
Mn spin, describes a transfer of spin from the Mn to the 

4(sk-,)^—n? +2nk-, 

electron system. The preractor - - -- - is zero 

for riki € {0,2}. For nk^ = 0 no transfer can occur be¬ 
cause there are no electrons that can exchange their spins 
with the Mn atoms; for nk^ = 2 the transfer vanishes due 
to Pauli blocking. 


The term proportional to (sk) x ((ski) x (S)) has the 
form of the relaxation term of a Landau-Lifshitz-Gilbert 
(LEG) equation and describes the tendency of a spin in a 
given effective magnetic held to align along the direction 
of the held. Unlike in the LEG equation, here, the pref¬ 
actor is determined by the parameters of the microscopic 
model and is not a phenomenological htting parameter. 

The last term in Eq. (10) resembles a relaxation term 
that would be expected in the EEG equation for the Mn 
magnetization (S). Here, it arises in the equation for 
the electron spin rehecting the conservation of total spin 
which is a feature also of the full QKT^®. However, there 
is a crucial difference between the last term in Eq. (10) 
and the EEG relaxation term for the Mn magnetization: 
while the cross-products in the EEG equation involve 
classical vectors, we are dealing here with vector oper¬ 
ators. Here, the expectation value has to be taken after 
constructing the cross product in a symmetrized form. 
The physical consequences of this difference become most 
obvious by rewriting the last term in Eq. (10) as follows: 

(S X (S X (sk))) + (((sk) X S) X S) ^ 

2 

- {(S^) - (5ll'))(sii) - + (5"'))(Si}), (12) 

where (sj}) and (si}) describe the electron spin of the 
states with k-vector k in the direction parallel and per¬ 
pendicular to the Mn spin vector (S) and S'!! = S • 

It is seen from Eq. (12) that even when the electron 
spin is aligned parallel to the Mn spin, a spin transfer can 
occur, and it was already noted in Ref. 16 that the cor¬ 
responding parallel spin transfer rate coincides with the 
result of Fermi’s golden rule. In contrast, the correspond¬ 
ing term in the standard EEG equation would be zero. 
This transfer is enabled because the factor (S^) — (S'!! ) 
is non-zero as quantum mechanically the maximal value 
of is while (S'^) = fi^S{S + 1) which reflects 

the uncertainty between the respective spin components. 
For classical vectors, as considered in the standard EEG 
equation, this factor would always be zero. Furthermore, 
in general the contribution in Eq. (12) is different for the 
parallel and perpendicular components of the electron 
spin. It is noteworthy that if the Mn spin had been rep¬ 
resented by a pseudo-spin i, this feature would be lost 
as then independent of the Mn spin conhguration we find 
{S"^) = I and ) = \ resulting in the same prefactors 
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for (s{[^) and (sj^^) in Eq. (12). 

In order to use Eq. (10) in practical calculations we 
have to know the values of the average Mn spin (S) and 
according to Eq. (12) the second moment ) which 
appear on the r.h.s. of Eq. (10). The average Mn spin 
can be calculated from the knowledge of the electron 
spin and the initial total spin by using the total spin 
conservation^®. Setting up an equation of motion for the 
second moment is cumbersome and not necessary for the 
cases that we shall discuss in this paper where it is as¬ 
sumed that the number of Mn ions by far exceeds the 
number of photo induced electrons {Nmu ^ Ne). In this 
case, the change of the average Mn spin as well as its sec¬ 
ond moment can be neglected and thus the second mo¬ 
ment essentially coincides with its initial thermal value. 
Furthermore, for nearly constant Mn magnetization, the 
equations of motion for electron states with different en¬ 
ergies ftwk are decoupled in the Markov limit due to the 
delta-distribution in Eq. (10) and the fact that re¬ 
mains constant which allows using the initial occupation 
for the evaluation of the frequency renormalization. 

The decoupling of the equations of motion in the 
Markov limit enables us to find analytical solutions for 
Eq. (10). To this end we split the electron spin into its 
components parallel and perpendicular to the Mn spin 
according to: 




■ COs(wMt) 


l(S)x(sk,(0))| 
((S)x(sk,(0)))x(S) 
|((S)x(sk,(0)))x(S)| 


(13) 


where lom accounts for the precession of the perpendic¬ 
ular component that results from Eq. (10). With this 
decomposition, Eq. (10) can be rewritten as: 


dt 


-7k,((52)-(5ll'))4.> 

= 7ki4iSki'S'- ^7ki((5'^) + {5'll^))si^^, 


with 


(14a) 

(14b) 


7ki = 


h^v 


k 


t^k), 


UlM 


Jsd'^Mn 

h 


S 



(15a) 

(15b) 


and S = |(S)|. Eq. (14a) is a Riccati differential equation 
with constant coefficients which can be solved analyti¬ 
cally. Its solutions can then be fed back into Eq. (14b) 
for the perpendicular electron spin. The explicit solu¬ 
tions are listed in appendix A. 

It is noteworthy that by a rescaling of the time axis 
according to r := 7kit all material parameters can be 
eliminated from Eqs. (14) for the moduli and 


Therefore, with this choice of time units and given initial 
conditions we obtain the same universal solution for all 
material parameters. Reinserting the solutions for 
and into Eq. (13) and choosing again l/7ki as the 
unit of time, we conclude that for given initial conditions 
the time trace of the electron spin (sk^) is affected by the 
material parameters only via the ratio ujm /7ki ■ 


III. NUMERICAL RESULTS 

The quantum kinetic equations of motion (7) have been 
solved numerically and compared with their Markov limit 
(10) for different initial conditions in a three dimensional 
bulk DMS. The initial electron distribution over the sin¬ 
gle particle energies Ek is taken to be Gaussian with its 
center at Ek=o and a standard deviation of cr = 3 meV 
while the initial magnitude of the Mn spin is set to 
(i.e., 20% of its maximal value). The material parame¬ 
ters used were the same as in Ref. 16 for Zno.gaMno.orSe 
with Aid = 12 meVnm® and mg = 0.21mo. 

First, we shall discuss results where at the beginning 
of the simulation the electron spins are assumed to be 
totally polarized in a direction with an angle of 45° with 
respect to the Mn magnetization vector. Displayed in 
Fig. 1 is the corresponding time evolution of the electron 
spin; part (a) shows the total electron spin, while in parts 

(b) and (c) the components parallel and perpendicular 
to the Mn magnetization are plotted, respectively. The 
full quantum kinetic results are plotted as solid red lines 
whereas curves derived from the analytical solutions of 
the Markov limit equations are depicted as dashed green 
lines. 

As seen from Fig. 1 (a), the dynamics predicted by 
the full theory is qualitatively different from the Marko¬ 
vian result. On a short time scale (for our parameters 
t < 5 ps), the electron spin decays much faster for the 
full solution than in the Markov limit. Subsequently, the 
quantum kinetic curve exhibits a non-monotonic time de¬ 
pendence and the electron spin eventually approaches a 
finite value. In contrast, in the Markov limit, we find a 
monotonic, almost exponential decay for all times. From 
the explicit analytical expression (cf. appendix A) it is 
seen that the long time limit of the electron spin in the 
Markov limit is zero. 

The origin of the non-monotonic behavior can be un¬ 
derstood by splitting the total electron spin into its com¬ 
ponents parallel [Fig. 1 (b)] and perpendicular [Fig. 1 

(c) ] to the Mn spin. Both spin components decrease al¬ 
most exponentially in the ML as well as in the full QKT. 
The time evolution of the perpendicular spin component 
essentially yields the same results for the full quantum ki¬ 
netic calculation and the Markov limit. In the full QKT, 
however, the parallel spin component changes its sign 
and converges to a finite negative value, whereas both 
spin components in the ML and the perpendicular spin 
component of the QKT drop to zero. When the parallel 
spin component in the full QKT crosses the zero line, its 
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time (ps) 


FIG. 1. (Color online) Time evolution of the total electron spin polarization (a) and its components parallel (b) and perpen¬ 
dicular (c) to the Mn spin assuming the electrons to be initially spin polarized along a direction at an angle of 45° relative 
to the Mn magnetization. The solid red line describes the spin dynamics according to the full quantum kinetic theory, the 
dashed green line shows its Markov limit (analytic solutions, cf. appendix A). Blue circles and purple squares correspond to 
approximate quantum kinetic calculations where only a subset of source terms for the correlations (as indicated in the key of 
the figure) has been accounted for. 


modulus has a minimum which leads to a minimum in 
the total spin. 

The obvious discrepancy between the different levels of 
theory with regard to the dynamics of the parallel spin 
component does not arise due to the assumption of a 
short memory in the ML. This can be seen from calcula¬ 
tions, where only the source terms , i.e., the terms 

used to derive the ML in the first place, have been taken 
into account but the finite memory expressed by the re¬ 
tardations in Eq. 8 are still kept [blue circles in Fig. 1]. 
The resulting curves almost coincide with the Marko¬ 
vian calculation. The main difference between the full 
QKT and the ML is due to the source term , 

which is demonstrated by simulations that incorporate 

only and [purple squares in Fig. 1]. 

The results of these calculations agree very well with the 
predictions of the full theory, suggesting that all other 
source terms are of minor importance, at least for the 
parameters used here. It should be noted, that especially 

the term \ like and b\XXf\ gives 

contributions to the reduced electron density matrices in 
the order of 0{Jg^) while the leading order contributions 
of the correlations are of 0(J^^). Thus, our results imply 
that a proper description of the coupled electron and Mn 
spin dynamics requires a treatment beyond perturbation 
theory. 

The effect of these higher order contributions on the 
dynamics is particularly dramatic in the case of initially 
unpolarized electron spins. Corresponding results are 
displayed in Fig. 2. Here, even the sign of the spin po- 



FIG. 2. (Golor online) Dynamics of the electron spin polar¬ 
ization for initially unpolarized electron spins. Line styles and 
symbols have the same meaning as in Fig. 1. 


larization is opposite for the QKT and ML calculations. 
Furthermore, also the predictions concerning the magni¬ 
tude of the spin polarization deviate significantly. 
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IV. INTERPRETATION OF THE SOURCE 
TERMS 


By the numerical analysis in the last section, we were 
able to trace back the difference between the full quan¬ 
tum kinetic theory and its Markov limit to a few selected 
source terms for the correlations in Eqs. (7). In this sec¬ 
tion, we shall give a physical interpretation to the indi¬ 
vidual source terms which will enable us to understand 
what determines their relative importance. 

First of all, ' is the most important source 

term, because it starts the correlation dynamics, i.e., 
without these sources the correlations would stay zero for 
all times. In the Markov limit, yields a Landau- 

Lifshitz-Gilbert-like damping term described in Eq. (12) 
and a spin transfer term proportional to the Mn spin 
(S). provides corrections for Pauli blocking to 

the transfer term and yields another LLG-like damping 
term, where the electron spin appears twice in the double 
cross product (cf. Eq. (10)). As seen above, the quan¬ 
tum kinetic contributions act similarly to their 

Markov limit counterparts. The dominant role of these 
terms is further emphasized by the fact that they are the 
leading terms in a perturbative treatment with respect 
to the exchange coupling constant Jsd- 

In order to understand the meaning of the 
terms, it is instructive to reformulate the equations of 
motion of the QKT by introducing new correlation func¬ 
tions according to: 


_\ ^ C/3 ^OC /-)^2'^2k2 

^ 01^1 ' / ^ ^nin2 *^1^2 ’ 

^1^2 

nin2 


(16) 


and eapj is the Levi-Civita symbol. We note in pass¬ 
ing that the residual sources contain a term re¬ 
sulting from which cannot be expressed by 

the summed correlations. Thus, Eqs. (17) are numer¬ 
ically advantageous only if is disregarded. 

The point here is that the two terms in Eq. (17b) orig¬ 
inating from and both involve the 

Levi-Givita symbol and can therefore be interpreted as 
describing precessions. This can be made more explicit, 
e.g., by introducing a vector with components a accord¬ 
ing to 

= (19) 

Then, the first of these terms, which stems from 
^hniki ’ 1*® written as a cross product: 

X (20) 

indicating a precession of the vector around the 

direction ujm of the Mn magnetization with the same 
frequency as the mean field precession of the electron 
spin. Likewise, the term originating from has 

a similar structure. It can also be written as a cross 
product 

X (21) 

where now the index j3 is associated with the components 

of a vector formed from the correlations according 

to 


(Qtj)/3 = QZ: 


cx.k.2 


( 22 ) 


which are summed over the electron band and Mn state 
indices. Here, we use the conventions a = 0... 3 with 
^hh ~ S'H'l /3 = 1...3. From Eq. (7c), we obtain 
the following equations of motion for the summed corre¬ 
lations: 

~ (17a) 




kX 

K,X 

(17b) 

where 


iJm =^nu„{S‘), 

(18a) 

k 

(18b) 

,ak 2 ^®® _ c/3 f; 

+ (18c) 

hh 

nin2 



i.e., now we are dealing with a precession around the 
direction of the electron spin. Thus, not only the 
average spins of the electrons and Mn atoms exhibit a 
precession dynamics, but also their correlations, which is 

represented in the equations of motion by the 
terms. 

Finally, the physical meaning of the source 

terms becomes clear by noting that their structure is 

analogous to the structure of the terms, where 

the products of electron and Mn density matrices are 
replaced by the corresponding unfactorized correlation 

functions. Thus, the sources provide the cor¬ 
related parts of the sources which represented 

a Landau-Lifshitz-Gilbert-like dynamics including Pauli 
blocking. 

Now that all source terms have been physically inter¬ 
preted, let us come back to the question of their relative 
importance in the case considered numerically in subsec¬ 
tion III. As already noted, the sources always 

play a pivotal role, since no correlations would build up 
without these terms. The importance of the remaining 
terms depends on the physical situation. Looking at the 


definition Eq. (7d-f) of the sources, it is seen that the 
terms , with X G {I, II, III}, comprise simi¬ 
lar factors as the corresponding contributions , 

except that the former contain an additional factor pro¬ 
portional to the electron density matrix . From this 

observation we can conclude that the sources 

X 1 

should be less important than the terms, if 

the electron density is moderate, as it is the case here. 
A criterion for being in the low density limit is partic¬ 
ularly easy to formulate for the terms, since 

Eq. (7e) implies that is negligible compared 

with if Nmtl ^ A"e which is fulfilled in our 

simulations. However, it is more challenging to give a 

condition for the negligibility of the term, as it 

strongly depends on the electron distribution in k-space. 

Finally, since the sources have the same 

structure as the term, except that the correla¬ 
tions place of the product 

they will be of minor importance if the relation 

gl2"2k2 

<C 1 is satisfied. The latter relation is expected 

‘iki ”■! 

to hold, when the conditions for the applicability of the 
correlation expansion are fulfilled. The numerical results 

shown in Fig. 1 indicate that the terms provide 

insignificant quantitative corrections which conhrms the 
consistency of the correlation expansion approach. 

The fact that a source contains correlations is, how¬ 
ever, not sufficient for concluding that it can be ne¬ 
glected compared with the terms, which do not 

involve correlations. In particular, the term 

was shown to qualitatively modify the spin dynamics 
(cf. Figs. 1 and 2). In view of our interpretation of 

the ^ term, this implies physically that account¬ 

ing for the precession of the correlations around the Mn 
magnetization is essential for a correct description of the 
spin dynamics. This also explains why previous studies 
in Refs. 16 and 17 reported a negligible contribution from 

the term, since there a situation was consid- 

ered, where the average Mn spin was initially set to zero 
which suppresses the precession. 

The features of the spin dynamics predicted in this ar¬ 
ticle manifest themselves in the time evolution of the spin 
polarization which is a quantity accessible experimen¬ 
tally, e.g., by time- and polarization-resolved photolu¬ 
minescence or Faraday-/Kerr-rotation measurements^^. 
Favorable for the observation of such effects should be 
experiments measuring the time dependence of the spin 
polarization as well as the its equilibrium value where the 
angle between the Mn magnetization and the initial elec¬ 
tron spin polarization induced by a circularly polarized 
laser beam is varied. For our purposes, bulk materials are 
preferable compared with, e.g. quantum wells, since for 


heterostructures, the anisotropy with respect to growth 
axis as well as structure inversion asymmetry can play a 
role^^ which would make it hard to separate the angular 
dependence predicted by our theory from anisotropy ef¬ 
fects. Furthermore, II-VI DMS should be better suited 
for the proposed experiment than III-V DMS, since they 
have the advantage of isoelectrical doping. In III-V ma¬ 
terials, the Bir-Aronov-Pikus interaction^^ between elec¬ 
tron and hole spins can dominate the spin dynamics®, 
while for II-VI DMS with sufficiently high Mn doping the 
s-d-exchange interaction is typically the most important 
spin relaxation mechanism^^. 


V. SUMMARY 

In this article, we have analyzed the spin dynamics of 
conduction band electrons in Mn doped bulk DMS in¬ 
duced by the s-d-exchange interaction. In contrast to 
our previous studies^®’^^, we now assume a non-zero Mn 
magnetization. This naturally leads to a distinction be¬ 
tween the electron spin dynamics of the components par¬ 
allel and perpendicular to the Mn spin which introduces 
an anisotropy in the spin relaxation. Starting from a 
microscopic quantum kinetic theory based on correla¬ 
tion expansion we have derived the Markov limit yield¬ 
ing equations similar to the widely used phenomenolog¬ 
ical Landau-Lifshitz-Gilbert equations. Our derivation 
yields microscopic expressions for the parameters in the 
Landau-Lifshitz-Gilbert equations and allows us to iden¬ 
tify some quantum corrections. The resulting rate equa¬ 
tions were solved analytically. 

Numerical simulations within the quantum kinetic the¬ 
ory revealed that, while the dynamics of the perpendic¬ 
ular electron spin component can be well described by 
the Markovian theory, the parallel component exhibits 
qualitative deviations between the full quantum kinetic 
and the corresponding Markovian results. The differ¬ 
ences between both levels of theory manifest themselves 
in a non-monotonic temporal behavior of the total spin 
in the quantum kinetic theory as opposed to an almost 
exponential monotonic decay predicted by the Markovian 
theory. Moreover, for certain excitation conditions, even 
the sign of the spin polarization differs between these 
levels of theory. 

A detailed analysis allowed us to assign a physical in¬ 
terpretation to all source terms for the correlations and 
to understand their relative importance found in our nu¬ 
merical studies. With the help of this analysis and our 
numerical results, the deviations between the full quan¬ 
tum kinetic theory and its Markov limit were traced back 
to the neglect of a precession dynamics of the correlations 
in the Markov theory. This precession is missing in the 
Markov limit not because of the assumption of a short 
memory but due to the perturbative treatment that is 
implicit in this approach. 



9 


ACKNOWLEDGMENTS 


The solution of Eq. (Al) is then given by: 


We acknowledge the support by the Deutsche 
Forschungsgemeinschaft through the Grant No. AX 
17/9-1. 


Appendix A: Analytical solutions of the Markov 
equations 


■Ski (f + 

where ip is determined by the initial value of s{[^. 

Eq. (14b) for the perpendicular spin component as¬ 
sumes the form: 


Eq. (14a) is a Riccati differential equation 



-?+/4i)ski> 


(A5) 


dt 


/4i -9si^+h, 


(Al) 


with / = 7 kiS', g = 7ki((-S'^) - (‘S'll^)) and h = 

7ki S ^ . For / = 0, which is the case if 5 = 0, 
the solution of Eq. (Al) is simply: 


where ^ = iyki + (^H^)^. Eq. (A5) is solved by 


■Ski (^) 


4i(0)e 


e 0 


=:/ 


(A6) 


4iW = (slii(0)-^)e-«‘ + 4 (A2) 

For / 0, the Riccati equation can be rewritten in terms 

of a linear differential equation with eigenvalues: 


For / = 0, / = 1 and the perpendicular spin component 
decreases exponentially. Inserting the solution for the 
parallel spin compontent from Eq. (A4) for non-zero / 
yields: 


\/2 



(A3) 


cosh (I) 
cosh + vt) 


(A7a) 
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